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This article is aimed at the investigation of some properties of the Weibull cumulative 
exposure model on multiple-step step-stress accelerated life test data. Although the model 
includes a probabilistic idea of Miner's rule in order to express the effect of cumulative 
damage in fatigue, our result shows that the application of only this is not sufficient to 
express degradation of specimens and the shape parameter must be larger than 1. For a 
\ random variable obeying the model, its average and standard deviation are investigated on 

a various sets of parameter values. In addition, a way of checking the validity of the model 
i is illustrated through an example of the maximum likelihood estimation on an actual data 

E~~ i ■ set, which is about time to breakdown of cross-linked polyethylene- insulated cables. 

1 Introduction 

S: 

In many industrial fields it is requested for lots of products to operate for a long period 
■ of time. Accompanied with that, it is very important to give reliability in relation to the 

lifetime of products. In such cases, however, life testing under a normal stress can lead to a 

lengthy procedure with expensive cost. As a means to cope with these problems, the study of 
[ accelerated life test (ALT) has been developed. The test makes it possible to quickly obtain 

information on the life distribution of products by inducing early failure with stronger stress 

than normal. 

One important way in ALT is step-stress accelerated life test (SSALT). There are mainly 



two types of SSALTs, a simple SSALT and a multiple-step SSALT. In the simple SSALT there 
is a single change of stress during the test. Miller and Nelson (1983) have shown optimum 
simple SSALT plans in an exponential cumulative exposure (CE) model. Xiong (1999) has 
studied an exponential CE model with a threshold parameter in the simple SSALT. Park 
and Yum (1998) have shown optimum modified simple SSALT plans in an exponential CE 
model, under the consideration that it is desirable to increase the stress at some finite rate. 
Lu and Rudy (2002) have dealt with the Weibull CE model with the inverse power law in 
the simple SSALT. 

On the other hand, in the multiple-step SSALT there are changes of stress more than 
once. Yeo and Tang (1999) have investigated a three-step SSALT in an exponential CE 
model. Khamis (1997) has proposed an exponential CE model with k explanatory variables 
and investigated it on three-step SSALT data. McSorley, Lu and Li (2002) have shown 
the properties of the maximum likelihood (ML) estimators of parameters in the Weibull CE 
model with a log-linear function of stress on three-step SSALT data. Nelson (1980, 1990) has 
proposed an important idea, which gives the basic CE model for life as a function of constant 
stress from SSALT data. This is a probabilistic analog of Miner's rule (Miner, 1945), which 
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is stated on a deterministic situation, and gives the basis of all models mentioned above. He 
also performed the ML estimation in the Weibull CE model with the inverse power law on 
multiple-step SSALT data concerning time to breakdown of an electrical insulation. Hirose 
(1996) has proposed a generalized Weibull CE model, which has a threshold parameter. 

As we have seen, there are many kinds of studies about SSALTs on the basis of the CE 
model and these provide significant understanding of the Weibull and exponential models 
and SSALTs. However, the validity of the models is not necessarily clear (Bagdanavicius, 
1978; Nelson, 1980, 1990). In this article we devote ourselves to considering the following 
questions: 

• Can the models really express degradation of products? 

• If so, what condition on the parameters is necessary for it? 

• When a random variable obeys the Weibull CE model with a threshold parameter, 
how do its average and standard deviation behave under a condition? 

In Section [2] we introduce the CE model with a threshold parameter, which is a gener- 
alization of the CE model provided by Nelson. After giving the Weibull CE model with a 
threshold parameter in Section [3], we analyze it in Section HI In Section [5] we give an ML 
estimation procedure. In Section [6] we illustrate an example of the ML estimation and a 
goodness of fit test on an actual data set and lastly give the conclusions. 

2 Cumulative exposure model 

We construct a generalized CE model with the help of the CE model proposed by Nelson 
(1980, 1990), whose model gives the distribution function of a random variable for failure 
time. Although the general model is obtained in a similar way to Nelson's, it differs in having 
a threshold parameter that decides whether a specimen is influenced by stress or not. 
The assumptions to obtain the CE model were given by Nelson as follows: 

i) The remaining life of specimens depends only on the current cumulative fraction ac- 



ii) If held at the current stress, survivors will fail according to the distribution function 
for that stress but starting at the previously accumulated fraction failed. 

Using a distribution function F of a non-negative random variable with an explanatory 
variable V and a threshold V t h, we construct the distribution function G of a random variable 
T for failure time in a sequential way. Denote by Vi a stress that a specimen is subjected to 
in an interval (t«-i, U] (i=l,2,. . . ). 

First of all, we define G by 



cumulated. 




(Vi > v th ) 
(Vi < V th ) 



(2. 1) 
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for to < t < t\. 

Next, for t\ < t < t 2 we define 




def J F(t-t 1 + S 1 ;V 2 ) (V 2 > Vth), 

' I F( Sl ;V 2 ) (V 2 <V th ). 



(2. 2) 



Here, according to Assumption ii), si is a positive value satisfying G(t\) = F(si\ V 2 ). 
Similarly, for t^i < t < tj we define 



where Sj_i is a positive value satisfying G(tj_i) = Vi). 

Also in Hirose (1996), a similar formulation was given, provided that the step stress at 
the present time is not lower than that at the past time, which means V m < Vth holds for 
any m < i when Vi < V t h- Actually, in his formulation 



We deal with a multiple-step SSALT under the condition that specimens were subjected to 
a normal level of stress and did not fail before the test. As we will see in the next section, 
this setting has a possibility of throwing light on new aspects concerning the Weibull and 
exponential CE models. In this section, first we introduce the multiple-step SSALT and 
second we give the Weibull CE model under the condition. 



During the multiple-step SSALT, specimens are subjected to successively higher levels of 
stress as follows. After a specimen was used at a normal level of stress, it is subjected to 
an initial level of stress for a predetermined time interval at the first stage in the test. If it 
does not fail, it is subjected to a higher level of stress for a predetermined time interval at 
the next stage. In analogy, it is repeatedly subjected to higher levels of stress until it fails. 
The other specimens are tested similarly. The pattern of stress levels and time intervals is 
the same for all specimens. 

3.2 Weibull CE Model with a threshold parameter 

We construct the Weibull CE model by combining the basic CE model in Section [2] and a 
Weibull law, in which the underlying distribution is the two parameter Weibull and the scale 





for < t < ti. Note that our formulation is more general. 



3 Model for SSALT 



3.1 Multiple-step SSALT 
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parameter is replaced with the function of an explanatory variable. In addition, we give the 
log likelihood function in the case when the step-stress data are given under the condition 
mentioned above. 

Let V be an explanatory variable and a function of it <p(V). When the scale parameter 
is replaced with <f)(V) in the Weibull distribution function, the distribution is given by 



We use this as F in Section [2j 

Denote by V s and T s a normal level of stress and the length of the time interval during 
that a specimen is used before the test, respectively. In addition, denote by Vi the stress 
that a specimen is subjected to at the (z — l)-st stage in the test, and let be the start 
time of the stage (i = 2, 3, 4, . . .). Since the level of stress becomes higher as the stage in the 
test advances, the relationship Vi < Vj holds when 2 < i < j. 

Before we consider the Weibull CE model under the condition mentioned in the first two 
lines of Section [31 as preliminaries, let us consider the model without assuming the condition 
holds. When we denote by ti the start time of the test and set to at t\ —T s , (12. lj) gives 







This and 




yield 




(V < Vth). 



Hence, (12. 2p gives 
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By repeating similar calculations, we can obtain the cumulative distribution function G in 
the (i — l)-st stage: 



G(t) = 1 -exp -e {t) 



U—i < t < ti 



(3. 1) 



where 



e(t) 



clef 



^(K) <j>(y k ) 



U-l — ti-2 

m-i) 



t - 1 



i-1 



tk — tk-l _|_ _|_ U-l — U-2 _j_ t — tj_i 



{v s > v th ), 

(V s < Vth) 



(3. 2) 



when 



V 2 < V 3 < ■ ■ ■ < Vfc_ x < Vth < V k < ■ ■ ■ < Vi-i < v. 



Since we are interested in the case of degradation of products under a normal level of stress, 
we assume V s > V t h in the sequel. 

Next, we seek our target, that is, the cumulative distribution function under the condition 
that a specimen was subjected to a normal level of stress and did not fail before the test. 
From the statements above, the function is as follows: 



G(t\t > tj. 



{G(t)-G(ti)}I {t>tl} (t) 

1 - G(h) 
{l - exp^ih) - e^(t)}} I {t>tl} (t), 



(3. 3) 



where 



ht>h}(t) 



def J 1 (t > t\ 



\ (t< tx). 

Thus, the log likelihood function In L under the condition is expressed by the following: if 
we denote by iV and I the sample size and the level of the stage at which a specimen fails, 
and use superscript (J) to show that a variable is related to the j-th specimen, 

TV N 

In^X^exp^tg^^ (3. 4) 

3=1 j=l 



where we express e(t) by e(t; T s ) in order to show clearly that each specimen has each T s 



4 Statistical properties 

We consider the statistical properties of the model under the condition mentioned in the 
previous section. First we state the role of the shape parameter /3 in the distribution function 
(13. 3p and second we investigate the relationship between statistical quantities and the values 
of parameters after we simplify the model without loss of generality. In the sequel we express 
G(t\t > ti) by G(t\t > ti;T s ) when it is necessary to show clearly the length of the time 
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interval during that a specimen is used before the test. In analogy, we express G(t\t > ti) 
by G(t\t > ti] V s ) when it is necessary to show clearly the normal stress that a specimen is 
subjected to before the test. Since we are interested in elapsed time from the start time of 
test, in the sequel we suppose that t\ is the base point in time. That is, we may consider t\ 
equal to 0. 

Depending on the magnitude of /3, the distribution function has a different aspect as 
follows. 

Lemma 4.1 Assume that T a < T b . Then, the following holds for t > tk-i such that Vk~i < 
V th < Vfc. 



i) If0<(3<l, 

ii) If 15 = 1, 
Hi) If f3 > 1, 



G{t\t>t 1 ;T a )>G{t\t>t 1 ;T b ). 
G{t\t>t 1 ;T a ) = G{t\t>t l ;T b ). 
G{t\t>t x] T a )<G{t\t>t 1 ;T b ). 



Proof. The substitutions of e(t\) = T s /<fi(y s ) and (13. 2h into (13. 3ft yield 



G(t\t > t X ]T s ) = 1 -exp 



j T s 1 j T s ^— v t m t m —i t t%—\ 



U(v a )j \m: 



+ E 



+ 



By differentiating this with respect to T s and arranging it, we find 



dG(t\t >h;T s ) (3 



Noting e(t) > e(ti), we can see 



e p -\t) - e^ih) exp e p {h)-£ p {t) 



i) if < < 1, 

G(t\t > ti, T s ) is a strictly decreasing function of T s since dG(t\t > t\\ T s )/dT s < 0, 

ii) if p = 1, 

G(t\t > t%\ T s ) does not depend on T s since dG(t\t > t\\ T s )/dT s = 0, 

iii) if p > 1, 

G(t\t > ti, T s ) is a strictly increasing function of T s since dG{t\t > tx\ T s )/dT s > 0. 
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This completes the proof. □ 
The statement i) in the lemma means that specimens become more durable as they are 
used longer before the test. This is clearly irrational. Thus, in this sense any value in (0, 1) is 
inadmissible for (3. The statement ii) deals with a situation when the underlying distribution 
is exponential. It indicates that the CE model inherits the memoryless property from the 
exponential distribution. The statement hi) expresses the most realistic situation, in which 
the durability of specimens decreases as the the duration of their use becomes longer before 
the test. 

In a similar fashion, we can obtain the following lemma. 

Lemma 4.2 Assume that (Vth <)Va < H- Then, the following holds for t > t^-i such that 
V k - X < V th < V k . 



i) //0</3< I, 

ii) If (3 = 1, 
Hi) If f3 > I, 



G(t\t>t 1 ;V a )>G{t\t>t 1 ;V b ). 
G(t\t>t 1 ;V a ) = G(t\t>t l ;V b ). 
G(t\t > ti,V a ) <G(t\t>t l ;V b ). 



From this lemma, we can know a similar fact to Lemma 14.11 Especially, note that the 
statement hi) expresses the most realistic situation, in which the durability of specimens 
decreases as the normal stress imposed before the test becomes higher. 

In the sequel we assume the inverse power law in for V > Vth' 

where K and n are positive parameters and Vth is a non-negative parameter. In addition, we 
assume that the length of the time interval and the breadth of upsurge of stress are constant 
in the test. That is, we set 

At ^ U - AV^V i+1 -Vi (i = 2,3,...) and V 2 = AV 

Let us simplify (13. 2p and seek the expectation and second moment of a random variable 
obeying (13. 3p . 

By using the above constants and rewriting (13. 2\i and ( 14. ip . we can obtain 
At ~ At At At t- ti_ x 

£{t) = W7) s + Wk) + "' + mT) + m) ^ ' (42) 

M C 1 "**)" At ((rn-l)AV-Vthr (;// ,., . , (4 . 3 ) 



^(K) k ' <p(y m ) k 
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where 



rp def T s 

s = At' 



^def 



A' 



(Aty s )«' 



f> def Vt/i 



and (A; — 2)AV^ < Vth < {k — 1)AV holds. Here, remark that A and Vth are parameters to 
be estimated while T s and AV^ are quantities to be prespecified in order to decide a concrete 
model. These expressions indicate that we can take At and V s as a unit of time and a unit 
of stress, respectively. Besides, we can suppose that < V t h < 1 when we deal with the case 
that V 8 > V th . 

By means of a similar procedure and the arrangement of expressions, we can obtain 
another e(t) for a different stress V/, say e'(t), in the following form: 



e'(t) 



At ~ At 

T s + — — + • • ■ + 



At 



+ 



At t-ti-x 



where 



<f>(V k ) <f)(V-i) <P{V t ) At 

At ((Vs/v s ) - v th ) n 



Note that only the first terms in the right-hand sides differ in the expressions of e(t) and 
s'(t). Thus, once we obtain the values of parameters, we can decide the distribution function 
in the case of another stress V' s (> V s ) by replacing only the first term in the right-hand side 
offlQI). 

Let us seek the expectation of a random variable T obeying (13. 3p . We first seek the 
following conditional expectation as preliminaries: for m > k, 



™ fU Q 

E[T\T < t m ] = Y, t^G(t\t > tx)t 



E 

i=k 



[tG(t\t > tx) 



U-1 



ti-i 



G(t\t > ti)t 



-t m exp[^(ti) - e\t m )\ + t k -! + E / 1 ex P[^(^i) - ^(t)]t. 



In the last line of this equation the relationship s(tx) = e(tk-x) is used, which holds by (13. 2D . 

When we denote by mo a positive integer such that l/<p(Vi) < 1 holds for any i > m , 
we can see that 



t m exp[— e (t m )] = (m — 1) At exp 

< (m — 1) At exp 
and the right-hand side converges to as m — > oo 



f A* t 

4>(Vi) J 



i=k 



f At ~ 2° At 
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Table 1: Parameter values 



k 


AV 


Vth 




n 


10 3 , 10 4 , 10 5 


0.39 


0, 0.5, 0.9 


2,3 


1, 2,3 



From the things above and E\T\ = lim E\T\T < t m ], we obtain 

m— ¥oo 

E[T] 



A/- Af Af . , 



Xt l=fc ^-i 

^E^r{-^) + ^-i)} (4-4) 



At At 

as the expectation in the case that At is used as a unit of time. Here, 



A(t) = expfe^ti) - ^(t)] | o °° I it + ^(t)} 1 



e u. 



The expression in the right-hand side of (14. 41) is useful for stable numerical calculations 
when T s takes a large value. 

In a similar fashion we obtain 

is f^y {-*&>+ «(*<-.)} 



At y /? v At ; 

as the second moment in the case that At is used as a unit of time. Here 
B t (t) d ^ expfe^tO-e^t)] 



(4. 5) 



e u. 



Using (14. 4p and (14. 5p . we can calculate the mean and standard deviation of T/At for the 
parameter values in Table 1. The results are shown on Figs 1 and 2. In these figures we 
show the difference in the pair of the values of Vth and n by the combination of the sort of 
line and the thickness or the color of line. That is, the solid, dash or dotted line means that 
Vth = 0, 0.5 or 0.9, respectively. On the other hand, the thick, normal or gray one means 
n — 1, 2 or 3, respectively. 
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Figure 1: The mean and standard deviation of T/At when /3 — 2 
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Figure 2: The mean and standard deviation of T/At when (3 = 3 
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5 Estimation procedure 



We state the way of seeking the ML estimates of the parameters in (13. 3p . (14. 2p and (14. 3p . 
Some techniques below help us to obtain the estimates numerically stably. 

We use a new parameter ( defined by K/K Q for a constant K instead of K. The reason 
is because only the parameter K possibly has an estimate that is much larger than those of 
the other parameters. 

By differentiating (13. 4p with respect to each parameter and arranging each equation, we 
can obtain the likelihood equations in the following simplified form: 

N N 

E^ + EM^^HO, 9e{P,n,C,V th }, (5.1) 

3=1 i 3=1 

where 

d * exp(-^(e;T^))-exp(-^(tp ) ;T^)), 

d = f f W) exp (-eP{t&; 2™)) + 5 9 (t®; W) t U) )) , 

S n (U,f s ) d ^ i^- 1 ^;^) ( t ((m - 1)A1/ - 4) n ln((m - 1)AK - 4)1 , 

^0 [m=k J 

5 c (t f ,f s ) ^ e^;T s ), ^(^^^^^(^^^^^((m-^Ay-^r- 1 . 

•"■0 m=k 

In the expressions above, note that e{t) in (14. 2p is expressed by e(t; T s ) as usual. We seek 
the zero of (15. ip by means of the damped Newton method (Bank and Rose, 1981) as follows. 

In this model, the calculation for ML estimates is so sensitive that, depending on a vector 
of initial guesses, a sequence of approximate vectors by the damped Newton iteration can 
converge to a vector of estimates in which the estimate of (3 is less 1 even if the value of 
the likelihood function is not a maximum value. According to our observation, when this 
phenomenon occurs, there is often a tendency that the estimate of V t h tends to 0. Thus, we 
adopt the strategy below. 

1. We seek the profile of (13. 4p with respect to V t h- That is, while changing the value of 
Vth from a value cto to another value ol\ in incremental steps, we seek the estimates of 
the other parameters in each step. 

2. Among the points on the part of the profile, we select the point at that the profile 
achieves its maximum, and then seek the ML estimates of all parameters simultaneously 
by using the point as a vector of initial values and performing the damped Newton 
iteration without fixing Vth- 

The derivatives of the expressions in the left-hand side of (15. ip are given in Appendix. 
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6 Example 



On the basis of the results obtained in the previous two sections, we show an example of the 
ML estimation in (13. 3p on an actual data set. The data set is part of the step-stress data on 
time to breakdown of cross-linked polyethylene-insulated cables in Hirose (1997). We chose 
only data whose insulation class is 22kV (in 3-phase). The reason is because data whose 
insulation class is 33kV in the literature indicate that their durability is higher than those 
for 22kV. Thus, we concluded that we can not mix both of them. Also, note Lemma 14.21 
Moreover, we did not incorporate into our sample data set the data coming from the cables 
that passed 26 years because their values are abnormal, compared with the others, and they 
push down the value of the likelihood function extremely. Finally, we perform a goodness of 
fit test by utilizing the Monte Carlo method. 

6.1 Data set 

We show the data set in Table 2. For each specimen j, the first column indicates the length 
of the time interval during that the specimen is used before the test, and the second column 
indicates elapsed time from the start time of test by the start time of the stage on which a 
specimen fails. The unit of time is ten minutes. The third column indicates the number of 
data in each row, which is denoted by Nd- The fourth and fifth columns indicate the average 
and standard deviation of data, respectively, in the case that an outrageous datum is not 
taken into account. 

6.2 Maximum likelihood estimation 

On the data set, we can see that AV = 0.39 because that V s = 22kV (in 3-phase) = 22/ v / 3kV 
(in single phase) and AV = 5kV (in single phase). As a constant for Kq, we set Kq = 10 4 
and used as a vector of initial guesses n, £, Vth) = (2.0,2,1,0.5). These were roughly 
guessed from the comparison between Fig. 1 or 2 and the averages or standard deviations in 
Table 2. In addition, for seeking the profile we set ao = 0.5, ot\ = 0.999 and the increment 
of V t h at 0.001 in each step. 

Ultimately, we can obtain the following ML estimates: 

P = 5.016812, n = 1.603875, ( = 0.548237, V th = 0.944054. (6. 1) 

Then, the value of the likelihood function (13. 4p is —244.4626. The mean and some statistical 
quantities of T and the scatter plot of data are given on Fig. 3. In the figure the solid line 
indicates the mean and the upper or lower dotted line indicates the mean plus or minus the 
standard deviation, respectively. Each dot indicates (t^}\ + t\ )/2 for each sample j. 
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Table 2: Step-stress data 



f(j) 

s 




N d 


Ave. 


SD 


157680 


54, 56, 59, 64 


4 


58.3 


4.3 


473040 


16 


1 


16 


* 


578160 


16, 16, 18, 19, 22, 24, 46, 48, 50 


9 


28.8 


14.7 


630720 


24, 24, 28, 30, 32 


5 


27.6 


3.6 


735840 


19, 20, 23, 35, 39, 42 


6 


30.0 


10.2 


788400 


12, 12, 12, 13, 14, 14, 14, 17, 17 


9 


13.9 


2.0 


840960 


12, 14, 14, 23, 24, 74+ 


6 


17.4 


5.6 


893520 


16, 18, 20, 20, 26, 35 


6 


22.5 


7.0 


946080 


16, 16, 16, 17, 18, 18, 20, 22, 26 


9 


18.8 


3.4 


998640 


8, 10, 11, 12, 12, 12, 13, 13, 14, 15, 22 


11 


12.9 


3.6 


1051200 


11, 12, 12 


3 


11.7 


0.6 


1156320 


11, 12, 12, 13, 14, 14 


6 


12.7 


1.2 



The marks f and * mean outrageous and incomputable, respectively. 
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2.6 25 50 75 100 
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Figure 3: Mean and scatter plot 



6.3 Goodness of fit test 

After the values of the parameters were estimated and a model was decided, we are often 
concerned with testing its validity. In the example, however, we can not perform a chi-square 
goodness of fit test on the grouped data because the number of samples is too small (Rao, 
2002, p. 396). Hence, in almost the same way as that Ross (2002, p. 206) adopted in such 
a situation, we test the hypothesis that the model is consistent with the data set when the 
parameters are of the values in (16. ip . 
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Simulation conditions 



We performed Monte Carlo simulation under the simulation conditions below. 

• Setting for generating simulated data 
We used as prespecified values 

AV = 0.39, K = 10 A , f s = 157680, 473040,..., 1156320, 

which are the same values as those for the data set in Table 2, and as the true values 
of the parameters 

P = 5.016812, n = 1.603875, C = 0.548237, V th = 0.944054, (6. 2) 

which come from (16. ip . 

For each value of T s , the number of simulated data is the same as that in Table 2. 
The outrageous datum in the table is, however, counted out. That is, for T s = 840960, 
the number of simulated data is 5 not 6. Hence, the total number of data is 74 in a 
sample set. 

• Sample sets 

Except sample sets where the ML estimates could not be obtained, 1000 sets of inde- 
pendent pseudo-random samples were considered. 

• Setting for estimation 

On the stage to seek the profile we set a = 0.85, cti = 0.999 and the increment of V t h 
at 0.001 in each step. 

Procedure for generating random sets 
From (GOD, ( EPD and fl4~3|) 

£(ti) ~m~ "IT" 

When we set u = f G(t\t > ii) for (oo >) t > this and (13. 3p give 



o-, 1//3 

e(t)= -ln(l-.)+(ii^f s ) 



(6. 3) 



Let us denote by q the expression in the right-hand side. Then, (14. 2p and (14. 3p yield 

g ((m-l)AV-V th r cq (l^%T fs< j2 ((rn-l)AV-V th ) n ^ (g 4) 

m=k -^oC -^oC m=k Ko( 

Note that u e (0, 1) and i > k since tk-i < t < oo. Consequently, the procedure for 
generating a pseudo-random set is as follows: 
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Table 3: Nonover lapping subintervals 



T 

J - s 


subinterval 


1 


2 


3 


4 


788400 


[0, 12] 


(12, 14] 


(14, oo) 




946080 


[0, 16] 


(16,18] 


(18, oo) 




998640 


[0,11] 


(11,12] 


(12,13] 


(13, oo) 



Table 4: Values of the test statistic 



T 

A s 


mi 


m 2 


m 3 


1714 


Pi 


P2 


P3 


P4 


T 


788400 


3 


4 


2 




0.122186 


0.067105 


0.810709 




26.22508 


946080 


3 


3 


3 




0.479264 


0.125719 


0.395017 




3.572318 


998640 


3 


3 


2 


3 


0.225639 


0.059853 


0.064528 


0.649980 


13.19004 



1. generate a uniform random number and seek q by (16. 3p . 

2. find i that satisfies (16. 4p . where the i is the stress level I in which failure occurs, 

3. calculate 

4. repeat the three steps above 74 times. 



Procedure for a goodness of fit test 

In Table 2 we choose some T s 's and calculate the value of the test statistic 

r = f^M (6 .5) 
£l N dpi 

for each of the T s 's. Here, stands for the number of data in the ith subinterval when the 
interval where failure time lies is divided into k nonoverlapping subintervals, and pi is the 
probability that failure time lies in the ith subinterval. 

In fact, we chose T s = 788400, 946080 and 998640 and for each of them divided the interval 
into subintervals shown in Table 3. The values of the test statistic and other variables in 
(16. 5p are, for example, as in Table 4 for the data in Table 2 and ( 16. 2p . 

In the simulation process we sought parameter estimates, calculated T on simulated data, 
and checked whether the value, say T sim , was at least as large as the value of T in Table 4. 
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Table 5: Successful number in satisfying T sim > T 



788400 946080 998640 simul. 



T ■ > T 

- 1 sim _ 



1 150 16 



Table 6: Biases and variances 



P n C 



Bias 
Variance 



0.201057 0.089020 0.833870 
0.496762 0.114639 7.444595 



-0.053946 
0.016145 



Simulation result 

The result is shown in Table 5. The second row indicates the number of the data sets on that 
^sim — T held. The last column indicates the number of the data sets on that the inequality 
held simultaneously in the three T s 's. From this, we can reject the hypothesis at any level 
a > 1 x 10 -3 because the p-value is less than 1 x 10~ 3 . 

The unsuccessful number was 158 in finding a vector of ML estimates. Both bias and 
variance of £ were relatively large as shown in Table 6. 

7 Conclusions 

Under the condition that specimens were subjected to a normal level of stress and did not 
fail before the test, we considered the two-parameter Weibull CE model with the threshold 
parameter in the multiple-step SSALT. This consideration revealed that the shape parameter 
/3 must be larger than 1 for the model to fit the realistic situation in which the durability of 
specimens decreases as they are used longer or with higher stress. 

After simplifying the model without loss of generality, for a various sets of parameter 
values we showed the average and standard deviation of failure time versus the duration T s 
of the specimen's use before the test. As we have seen in Section [61 we can utilize them in 
calculations for ML estimates. 

In the section we performed a goodness of fit test by means of Monte Carlo simulation. 
In general it is not clear whether the inverse power law holds for every stress appearing in the 
step-stress test, and it is not clear even whether the basic idea of the CE model is available. 
Including these things, the example we showed gives a way of checking the validity of the 
CE model. 
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Appendix: Derivatives needed for the ML estimation 



For #1, 9 2 G {(3,n,(,V t h}, the derivatives of the expressions in the left-hand side of (15. II) 
are given as follows: 



_d_ 

W 2 



where 



clef 



A 



(j) 



1 gAg 

dj 86 2 



N N 

E^ + EM*?;^) 

j=i a i i=i 

iV 

E 

iV 



+ 



(9^2 / dj ^ (96*2 9t9 2 



a2' } A# 1 OA? 



dj dj dj 9t9 



1 (0 = /3) 

f (* = n) 

-f (0 = 

-f (0 = ft 



dcf 



A' ' ;, 'exp(^(tp ) ;^))+^(tg 1 ;^))) 



t/ij 



(9 



5 



(1 - Ej 



In the above, note that the following equation does not necessarily hold because (15. ip is a 
simplified equation not the original likelihood equation: 



d lAA^ v 



8 \JL\i j) v 



#02 I £l rf i J=i 



Each derivative of 5g (9 E {(3,n(,V t h\) is given in the following. The arguments in some 
derivatives are omitted as far as it does not cause a confusion. 



W = 6 " ln£ ~ 

d6 p ' 



C f lne - - j <J C , 



dV t 



tit 



C-c r lne — - | 5t> , = (Line, 



d_ 

dn 



C n - e^- 2 (t t ;f s ) |i ^ ((m - 1)AV - Vk) n ln((m - 1)AV - V^)| 
+ e^fe; T.) (i- ((m - l)AV - V th ) n (ln((m - 1)AV - V th )f) , 



■£r s »&> f >) = ( c v th + ^ ) e/,_2 (**; ^) \y o E (("» - ^ - 1 



x i E((™ - X ) A ^ - ^) n M(™ - !) A ^ - ^) - 



ne^fe; f a ) \ -i- E ((rn - 1)AV - t4) n_1 ln((m - 1) AV - 14,) 

l-^O m =k 



C n - ±) e^\U- T s )(i^ ((m - 1) AV - V^f ln((m - 1) AV - t^) 



X 



E ((m - 1) AF - VtuY' 1 



rrt=fc ) 

+ e^ft; f s ) (i- £ ((™ - 1)AT/ - y,)- 1 ln((m - 1)AV - V^)) , 



if, 



m=fc 

= (c^ + ^ 2 (t,;T s ) {i E ((m - 1)AV - v^r 1 



l m =fc 



1 



- (n - l^fe; — E (("» " W " ^) n " 

I ^0 m =k 
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